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Abstract. Magnetoacoustic tomography with magnetic induction (MAT-MI) is a 
coupled-physics medical imaging modality for determining conductivity distribution 
in biological tissue. The capability of MAT-MI to provide high resolution images has 
been demonstrated experimentally. MAT-MI involves two steps. The first step is a 
well-posed inverse source problem for acoustic wave equation, which has been well 
studied in the literature. This paper concerns mathematical analysis of the second 
step, a quantitative reconstruction of the conductivity from knowledge of the internal 
data recovered in the first step, using techniques such as time reversal. The problem 
is modeled by a system derived from Maxwell’s equations. We show that a single 
internal data determines the conductivity. A global Lipschitz type stability estimate 
is obtained. A numerical approach for recovering the conductivity is proposed and 
results from computational experiments are presented. 
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1. Introduction 


Electrical conductivity of the biological tissues can provide important information for 
clinical and research purposes. Conductivity imaging has been a subject of research for 
decades and the literature is vast. 

Magnetoacoustic tomography with magnetic induction (MAT-MI) is a new 
noninvasive modality for imaging electrical conductivity distribution of biological tissue 
[mini US]. In the experiments, the biological tissue is placed in a static magnetic held. 
A pulsed magnetic held is applied and induces an eddy current inside the conductive 
tissue. Consequently, the Lorentz force, the force acting on currents in the static 
magnetic held, causes vibrations and the tissue emits ultrasound waves. The ultrasonic 
signals are measured around the object. MAT-MI belongs to the class of coupled-physics 
imaging method which is often refered to as ‘hybrid imaging’. For a review on hybrid 
imaging methods that recover electrical conductivity distribution, we refer to HU. 

Hybrid imaging typically involves two inverse problems. In MAT-MI the two steps 
are decoupled. The hrst step involves an inverse source problem for the acoustic 
wave equation. This problem has been studied extensively in many works including 
p El El EHl [la E]. The second step, the focus of this work, is to reconstruct the 
spatially varying electrical conductivity from knowledge of the acoustic source. 

In the MAT-MI experiment, the object to be imaged is placed in a constant static 
magnetic background held Bq = (0, 0,1). A pulsed magnetic stimulation is introduced. 
The pulsed held is of the form BiM(f), where the vector held Bi is a constant and u{t) 
is the time variation. The magnetic permeability of biological tissue is approximately 
equal to that of a vacuum. Therefore the tissue does not have any noticeable ehect 
on the magnetic held itself. As a result, the time-dependence of the electromagnetic 
helds is u(t) and we need only to consider their spatial dependence. Because the electric 
held will depend on conductivity a, we write it as Eg-. Let denote the domain to be 
imaged. Then it can be shown that the electric held satishes 



in fl, 
in fl, 
on dn 


( 1 . 1 ) 


The hrst step in the MAT-MI inverse problem is to recover the acoustic source in the 
scalar wave equation from observed data at a set of locations. The acoustic source is 
related to the electromagnetic held; knowledge of the acoustic source in this model is 
equivalent to knowing the quantity V ■ (crEo- x Bq) throughout 

In this paper, we focus on the second step of MAT-MI, i.e., reconstruction of the 
conductivity a from the internal data given by V ■ (cxEo- x Bq). Our main result is that, 
if the conductivity is a priori known near the boundary, then it can be uniquely and 
stably reconstructed from one internal data. More precisely, the main result of this work 
reads as follows. 


Analysis of the Magneto-acoustic Tomography with Magnetic Induction (MAT-MI) 3 


Theorem 1.1 Denote the forward map, the map from conductivity to acoustic source, 
as F{a) := V ■ x Bq) 
support of cTi 


Suppose that ai and 02 satisfy Assumption 2.1 and the 


02 is away from the boundary of Vt at a distance greater than some 
constant ro > 0. Then, there exists a constant K, which only depends on vq, \, A and 
12, such that, if 


then the ineguality 

ll^"! ~ <^21^2(0) < 4 ||F((Ti) — F{a2)\\L'^(n), 

holds true. 


( 1 . 2 ) 


(1.3) 


During the completion of this work, we discovered a recent paper by Ammari, 
Boulier and Millien [1]. Their work also focused on the conductivity reconstruction 
aspect of MAT-MI. What is different is that the authors chose to reconstruct hrst the 
current density in the medium. They propose methods to solve for conductivity from 
current density. In our approach, we directly deal with the relationship between the 
acoustic source and the electromagnetic held, and propose a method that hnds the 
conductivity from the acoustic source. 

The rest of the paper is organized as follows. Section 2 introduces the notation 
used and basic results needed. In Section we study the mathematical model of the 
second step of MAT-MI and the linearized version of this problem. Section]^ is devoted 
to addressing the uniqueness and stability estimate of both linearized and nonlinear 
problems. In Section 5, we propose an numerical method for solving the inverse problem 
and present some results from computational experiments. A hnal section discusses our 
hndings. 


2. Notations and preliminaries 


We begin by introducing the notations for the the mathematical analysis. Throughout 
this paper, the standard notations for continuous differentiable function spaces and 
Sobolev spaces are used. Let 12 be a bounded domain in with Lipschitz boundary 
dfl. A typical point x = {xi,X 2 ,x^) G denotes the spatial variable. We use the 
notation (7°° (12) for inhnitely differentiable functions on 12 and (7^(12) is a subset of 
which contains the functions with compact support. We use (■, •) to denote the 
inner product in the Hilbert space L^(f2). For p > 1, we denote by IF^’^(f2) the L^-based 
Sobolev spaces on 12 with the usual norm. 




n=l 


du 

dXn 


Lp{n) 


In the case p = 2, we use the notation H^{Q) = IF^’^(f2), which is a Hilbert space. The 
Sobolev space is dehned as the closure of (7“ (12) in The dual space of 

II^{Q) is denoted by If there is no danger of confusion, we omit the domain 12 
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and abbreviate with L^, H^, and In the following, we do not distinguish 

in the notation for inner product, function spaces and the corresponding norms between 
scalar- and vector-valued functions. 


Assumption 2.1 Let a be a positive function belongs to W^’°° and assume that 

a(x) > A, Vx G 12. (2.1) 

and 


ll^'llwpoo < A 
for some constants A,A > 0. 


We start with stating several useful results on the elliptic partial differential 
equations with Neumann boundary condition. 


Definition 2.2 We say that u G is a weak solution of the Neumann boundary value 
problem, 


if 



V ■ (ctVm) 

= -V ■ E, 

in 12, 

(aVu -\-E) ■ v 

= 0, 

on dfl, 

aVu ■ V(p dx = 

— / E ■ V(pdx, 
Jn 

V(p G 


( 2 . 2 ) 


We need the following regularity result and standard energy estimate of the 
gradient. 


Proposition 2.1 Suppose that a satisfies Assumption 2^. For field E G the 
Neumann problem (2.2) has a solution u G H^. The solution u is unigue up to an 
additive constant and satisfies the estimate, 

\\Vu\\l2<X-^\\E\ 


|L2. 


(2.3) 


Proof The proof of the existence and uniqueness up to an additive constant is a 
standard result by the Lax-Milgram Theorem. We refer the readers to [16]. In the 
following, we prove the gradient estimate (2.3). 


It follows from the ellipticity condition (2.1) that 
^ [ o'\Vu\‘^dx. 


Taking the test function (f in Definition 2.2 to be the solution u, we have that 


ctVm ■ Vu dx = 


E ■ Vudx. 


Consequently, applying the Cauchy-Schwarz inequality, we obtain that 


A||Vni|i2 < 


E ■ Vu dx 


< l|V“lln2||E||L2, 


and (2.3) follows. 


□ 
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3. Analysis of the forward problem 

3.1. The forward problem 


The second step of MAT-MI is modeled by (1.1), where u is the unit outer normal 
vector of dVt and Bi = (0,0,1) is a constant vector. The data for this inverse 
problem is the acoustic source recovered from the first step, namely, V ■ (crEo- x Bo) 
with Bq = (0,0,1). The inverse problem of the second step of MAT-MI consists of 
reconstruction of conductivity a from knowledge of V ■ (cxEo- x Bo). 

We refer the readers to p] for the regularity results of the Maxwell’s equations. In 


Proposition 3.1, we show some regularity results of our reduced system (1.1). 
Definition 3.1 We say that E^. E is a weak solution of the (1.1) if 


and 


[ E,, ■ (V X $) dx = [ $ ■ Bi dx, G 

Jn Jn 

/ crEo-■ V(p dx = 0, \/ipEH^. 

Jn 


(3.1) 


We dehne the forward problem as 

F : ^ 

cr !->• V ■ (ctEo- X Bo). 

Next, we introduce a proposition on the existence, uniqueness and uniform L^- 
boundedness of the electrical field E(j. This implies that forward operator F is well- 
defined. 


Proposition 3.1 Let a satisfy Assumption 2.1. Then the system (1.1) is uniguely 


solvable and there exists a constant C\ depending on X, A and H, such that 

||Eo-||i,2 < Cl. 


Proof This proposition will be derived as a consequence of Proposition 2.1 Let us hrst 


reduce the system (1.1) to a Neumann boundary problem. Let E = W—y,x,0)- We can 


readily check that V x E = Bi. Hence V x (Eo- —E) = 0 and we can write E^- = E-|-Vm. 


Substituting this into (1.1), we have that u solves the Neumann boundary problem. 


(■jE). 


in fl, 
on dVt. 


(3,2) 


V ■ (ctVm) = —V 
{aVu + (jE) ■ v =0, 

The existence of u and uniqueness of Vm follows from Proposition |2.1[ For the uniqueness 
of Eg-, we consider the equations 


V ■ (aVn) = 0, in H, 

aVn ■ V =0, on dVL. 


(3,3) 
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If both El and E 2 are solutions to the system (1.1), then we have that Ei — E 2 = Vn 


and V solves the equations (3.3). By Proposition 2.1, the only solutions to (3.3) are 


constants. Hence Vv vanishes and E^- is unique. 

What remains is to show the boundedness of E^. Applying Proposition 2H to 
u, we have that 

||Vm||l 2 < A"^||crE||i2. 

Hence, 

l|Eo-||i;,2 = ||E + VM||i2 < (A/A + 1)||E||i^2. 

Note that we can choose E = ^{—y -\-a,x-\-b,0) and repeat the above argument for any 
constants a and b. It follows that, 

||Eo-||l2 < Cl, 

where 


C'l = 7 ^(A/A + l)inf ||(-?/ + a,a: + 6,0)||L2, 


only depends on A, A and H. 


□ 


3.2. Linearizaton of the forward map 

Recall that the distribution of the electric field E^- depends nonlinearly on the 
conductivity a and V ■ (ctEo- x Bo ) is nonlinear with respect to a. It is natural to start by 
linearizing the relationship between conductivity and data. In this section, we introduce 
the linearized of the inverse problem. We first examine the Frechet differentiability of 
the forward operator F. Then, some useful properties of the Frechet derivative at a, 
DF^, are presented. 

To introduce the Frechet derivative, we consider the following Neumann boundary 
problem. 


V-(crV(ph) =-V-(hE^), 

+ hE^) ■ z/ =0, 

where h G is the increment to the conductivity. 


in n, 
on dVt, 


(3.4) 


Theorem 3.2 For a satisfying Assumption 2.1, the forward operator F, defined in 


is bounded and Frechet differentiable at a. Its Frechet derivative at a, DF^, is 
DF^{h) = V ■ {{aV^ph +hE^) xBo), (3.5) 


given by 


where ifh solves (S.f), and satisfies 

\\DF,{h)\\L2 < C'2||h||iyi,oc, Vh G 

for some constant C 2 depends on A, A and H. 


(3.6) 
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Proof We first prove the boundedness of F. We can write 
F(cr) = V ■ (crE,^ X Bq) 

= crV ■ (E^ X Bo) + Va ■ (E<^ x Bq) = a + Va ■ (E<^ x Bq). 


It follows, by boundedness of a and Proposition |3.1[ that 

\\F{a)\\L2 < \\a\\L2 + \\Va\\L^\\E4L^ < + Ci)\\a\\wi,^, 

where Ci is the same constant as in Proposition |3.1[ 

Next, we show the Frechet differentiability of F at a. Consider the data 

F[a + h) = V ■ ((<7 + h)Eo-_|_/j x Bq) 


for some h G such that a + h also satisfies Assumption 2.1, where is the 

solution to (1.1) with a replaced hy a + h. Note that 

V X - E,,) = 0. 

Hence we can write Ecr+h — E^ = Vn. Substituting this into the equations for E 0-+/1 and 
Eo-, we obtain that u solves 

V-(crVM) =-V ■ (hE^+/j), in 11, 

(aVn + KEfjj^h) ■ = 0, on dVt. 


(3.7) 


Applying Proposition 2.1 to u, we have 
\\Vu\\l 2 < X-^\\hE^+h\\L^. 


Let V = u — iph, where iph solves (3.4). Then, v solves 

V-(cTVn) =-V-(Wm), 

(aVn + hVn) ■ z/ =0, 


in n, 
on dH. 


Applying Proposition |2.1| to v, we have 
llVnlU^ < X-^\\hVu\\L2. 

To estimate the remainder terms, we write 

F{a + h) - F{a) - V ■ {{aV(ph + hE„) x Bq) 

= V ■ ((cr(Eo-+h — Efj — V(ph) + h{Erj+h — Ecr)) X Bq 
= V ■ {{aVv + hVu) X Bq) 

= Vcr ■ (Vn X Bq) + Vh ■ (Vu x Bq). 


(3.8) 


(3.9) 


(3.10) 


Therefore, by (3.8), (3.10) and Proposition 3.1, we have 
||F(cr + h)- F{a) - V ■ {{aViph + hE,) 
= II V(J ■ (Vn X Bq) + Vh ■ (Vn x Bo)||l 2 
< ||Vcr|Uoo||Vn|U2 + ||Vh|Uoo||VniU2 


X B 


oJI|l2 


< CiAA 


-2 


L°° ■ 


i^co+CiA-lV/ilUo 

We can readily check the linearity of the operator maps h to V ■ + hEg.) x Bo). 

This complete the proof of Frechet differentiability of F at a. 
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What remains is to show that the formal Frechet derivative DF^j is a bounded linear 
operator. Note that 

DF„{h) = V ■ {{crViph + hEa) x Bq) 

= CtV ■ X Bg) + V(T ■ (Wiph X Bg) + hV ■ (E^ X Bg) + Vh ■ (Eg- X Bg) 

= /l -|- VCT ■ Ph X Bg) F V/i ■ (E^ X Bg). 


By applying Proposition 2M to ph and Proposition ^ to E^r, we conclude that 
\\DF^{h)\\L2 

= ||h + Vcr ■ {Vph X Bg) + Vh ■ (E^ x Bg)||i2 

— II^IU^ + A|| V</9/i||i2 + II V/i||loo ||Eo-||l2 

< \n\^/^h\\L^ + Ax-^\\hEjL^ + CiiivhiUoo 

< (|fi|i/2 + C'i(AA-i + l))||h||vci.oo. 


□ 


4. Uniqueness and stability 

In the following theorem, we obtain a Lipschitz type stability estimate for the inverse 
problem under certain conditions on the conductivity. The uniqueness of the inverse 
problem follows. 


Theorem 4.1 Suppose that a satisfy Assumption 2^. If a only depends on the third 
component of the spatial variable, i.e, a(x) = cr(x 3 ), then the ineguality 


\DF„{h)\\L2i^^ > -||h||i2(f^) 


holds true for any h G hFo ’“(r2). 

Proof Note that 

V ■ {Vph X Bg) = 0, 

for any function ph and that 

dcT 

Va X Bg = (0, 0, X (0, 0,1) = 0. 

0 x 3 

Hence, 

V ■ {aVph X Bg) = aV ■ {Vph x Bg) + Va ■ {Vph x Bg) = 0. 

Therefore, 

DF„{h) = V ■ {{aVph + hE„) x Bg) = V ■ {hE„ x Bg). 
Multiplying the both sides by h and integrating over hi, we obtain that 

/ hDFa{h) dx= hV ■ (hEo- x Bg) dx. 

Jn Jn 


(4.1) 
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By using the integration by parts twice, we have 


hDF„{h) dx= hV ■ (hEg. x Bq) da; 


= — (hEo- X Bo) ■ Vhda; 


in 


[ (E. xBo)-V(h2)da; 
Jn 


= - / W-(E. xBo)da; 

^ Jn 

The last identity above follows by noting 

V ■ (Eo- X Bq) = V X Eg- ■ Bq = Bi ■ Bq = 1. 


Then, by applying Cauchy-Schwarz inequality to hDF„{h) da;, we obtain (4.1). □ 


The same technique can be used to provide a general stability estimate for the 
nonlinear inverse problem. Note that, in the following theorem and corollary, no 
smallness constraint on the difference of conductivities is needed. 


Theorem 4.2 

and 


Suppose that ai and a 2 satisfy Assumption 2.1 //cxi — (T 2 G iyo’“(r2) 


Veil X Vcr2 ■ Bq = 0, 
then the ineguality 

||o-i - cr2||L2(o) < 2 ||F((Ti) - F{a2)\\L^(n), 

holds true. 


(4.2) 

(4.3) 


(4.4) 


Proof Assume that Ei and E2 solve (1.1) with a replaced by Ui and (J2, respectively. 
Let us multiply F{ai) — F{a 2 ) by di — 02 and integrate over hi to obtain 

[ [cTi - o- 2 )(E(cri) - F{a 2 )) da; 

Jn 

= / (di - d 2 )(V ■ (diEi - d 2 E 2 ) X Bo) da; 

Jn 

= [ (ai- d 2 )(V ■ ((di - d 2 )Ei X Bo) V ■ (d 2 (Ei - E 2 ) x Bq)) dx 
Jn 

= 0 Iki - cr 2 ||i 2 (m + {(Ti- d 2 )V ■ (d 2 (Ei - E 2 ) X Bo) da;. 

Jn 

In the above inequalities, the last step follows by the similar argument as in the proof 
of Theorem 14.11 

Next, we estimate 

[ [cTi - 0 ' 2 )V ■ (d 2 (Ei - E 2 ) X Bo) da;. 
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Recall that V x (Ei —E2) = 0 . Hence, we can write Ei —E2 = Vm. Applying integration 
by parts twice, we obtain that 

[ (cTi - (T2) V ■ (cr2(Ei - E2) X Bo) dx 
Jn 

= - a2V(ai - (T2) ■ (Vm X Bq) dx 

Jn 

= / cr 2 (V(cri - (T2) X Bo) ■ Vndx 
Jn 

= - / V ■ (ct2(V(cti - CT2) X Bo))nda: 

Jn 

= - (0-2V ■ (V(cri - 0-2) X Bo) + Va 2 X V(cri - <72) ■ Bo)Mda; 

Jn 

= 0 . 


(4.5) 


Here we nse the equalities (4.2), 


and 


V ■ (V(ai - 0 - 2 ) X Bo) = 0, 


Vcr2 X Vct2 = 0. 


Combining (4.4) and (4.5), we discover 

/ - (T2){.F{ai) - F{a2)) dx = ^||ai - ^211^2(^2)- 


The stability estimate (4.3) follows by applying the Cauchy-Schwarz inequality to the 
left-hand side of the above equality. □ 


In the following corollary, we list some simple cases, in which, the criteria (4.2) is 
easy to check. 


Corollary 4.1 Suppose that cxi and 02 satisfy Assumption 
and satisfy any one of the following three conditions: 


2.1 


Ifa,-a2eW^’^{n) 


(i) (Ti only depends on the third component of the spatial variable x^; 

(a) There exists a real number t such that tai + {1 — t)a 2 only depends on x^; 

(Hi) There exist a positive integer N and real numbers On, n = 1,2...N such that 
<^2 + an{cri - CT 2 )” only depends on x^; 

then the stability estimate (f.tJ) holds true. 

Proof We can readily see that (i) and {ii) are simple cases of {Hi). It suffices to show 


that (4.2) is satished and apply Theorem 4.2 


From {Hi), we know that 


N 


V (72 + ^ 0„(<Ti - (T 2 )” X Bj = 0. 


(4.6) 


n=l 
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In light of (4.6) and the facts that 

Va2 X Va2 = 0, 

we have the following equalities, 

Veil X V(T 2 ■ Bq 

= V(cri - (T 2 ) X Vct 2 ■ Bo 


N 


= V(ai - 0 - 2 ) X V 1^0-2 - 1^0-2 + - 0 - 2 )' 

= V(cri - da) X j ^ a„n((Ti - era)"" V(cti - era) 

\n=l / 

' N 

ann{ai - a2)"'~^V{ai - aa) x V((Ji - exa) 


Bn 


Br 


Br 


vn=l 


= 0 . 


The proof is completed by applying Theorem 4.2 


□ 


Roughly speaking, in Theorem |4.2| , we prove that, if the structure of two 
conductivities satishes the condition (4.2), the inverse problem bears a Lipschitz stability 
estimate. We propose next to remove this structure condition. In Theorem ]_T, we show 
that, if one conductivity varies less dramatically, the Lipschitz type stability estimates 
also holds true. 

Proof of Theorem 11.11 The proof differs from the one of Theorem |4.2| in the treatment 
of the last term in (4.4), 


I — / (cxi - 0 - 2 ) V ■ (a 2 (Ei - Ea) x Bq) dx. 
Jn 


We continue from ( 4.4| ). First, we estimate the electric held difference. Note that Ei —Ea 
is curl-free and we set 

Wu = El — Ea. 

Then, u satishes the equation 

V-(criVM) =-V ■ ((cxi - (T2)E2), 

Vu ■ z/ =0, 


in 12, 
on dH. 


(4.7) 


Applying Proposition |2.1| to u, we obtain that 
||Vm||l 2 < A"^||(ai - a2)E2||L2. 

From the standard estimate of elliptic equations [3, Chapter 9] and the Sobolev 
Embedding Theorem, we know that Ea is bounded and 

||E2||l°° < Ci 

where C only depends on tq. A, A and 12. Thus, we conclude that 

||Ei — E 2 IIL 2 < C||(cri — cr2)||i;^2. 
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Now, with the choice of K such that KC < 1/4, we estimate |/| as follows: 




[ (c^i - cr 2 )V ■ (a 2 (Ei - E 2 ) x Bq) dx 
Jn 


' (cTi - 0-2) Vcr2 ■ ((El - E2) X Bo) dx 

n 

< ||Vcr2||L°° llc^l ~ 0'2 ||l2 ||Ei — E 2 IIL 2 

< C'||Vcr2||L- ||(Ti - 0-211^2 


Substituting (4.8) into (4.4), we discover that 


(ffi -CT2)(f((Tl) - F((T2))dl > -||(T, -CT2l|!a(n)- 


(4.8) 


The desired estimate (1.3) follows by applying the Cauchy-Schwarz inequality to the 
left-hand side. □ 


5. An iterative reconstruction scheme 

One possible approach to solving the inverse problem is to formulate it as a least- 
squares problem. One can then apply a gradient-based method to solve the least- 
squares problem. Such a method will require knowledge of the Frechet derivative of 
the forward map which we studied in Section Convergence analysis of this type of 
reconstruction approach is available in HE]. Results in these references, together with 
our analysis of of DF in Sections and can be used to provide a convergence analysis 
for the iterative reconstruction of MAT-MI using steepest descent method. The main 
challenge of the least-squares approach lies in the difficulty to accurately evaluate DF 
and its adjoint where numerical differentiations are involved. We temporarily abandon 
the least-squares approach in favor of one that is based on a hxed point method. This 
approach is described next. 

5.1. Formulation 

In view of the structure of this inverse problem, we propose a novel iterative scheme, in 
which, the forward map and its derivative are not required. The desired conductivity 
is updated by solving a stationary advection-diffusion equation. Let al denote the 
unknown conductivity to be reconstructed, El be the corresponding electric held and g 
be the internal data obtained in the hrst step of MAT-MI. The internal data is related 
to the conductivity and the held through 

(7 = V-(alEtxBo). 

The algorithm proceeds as follows: 

(SO) Select an initial conductivity Uo and set A: = 0; 
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(SI) Calculate the associated electric field by solving the boundary value problem 

= Bi, 


(5.1) 


V X Efc = Bi, in fl, 

V ■ (cTfcEfc) = 0, in Q, 

cTfcEfc • z/ =0, on dQ. 

(52) Calculate the updated conductivity by solving the stationary advection-diffusion 
equation: 

V ■ (cjfc+iEfc X Bo) = g, in 

(jfc+i = (To, on (diZ. 

(53) Set fc = /c + 1 and go to (SI). 

Convergence test can debased on ||(Tfc—crfc_i|| or based on data misfit H^f —V-((Jfc+iEfc+i x 

Bo)l|. 

5.2. Convergence analysis 

The main advantage of this scheme is two-fold: First, the update of the conductivity is 
calculated directly using the the measured data and the simulated electric field. Hence, 
fewer numerical differentiations are involved when compared to the gradient-based least- 
squares minimization. Second, the convergence analysis can be carried out using an idea 


similar to the one in the proof of Theorem ]_T A global convergence result and a linear 
convergence rate are established the following theorem. 

Theorem 5.1 Suppose that the true conductivity satisfies Assumption 2^ and 

llVa^llioo < 2K, (5.3) 

where the constant K is the same as in Theorem [73 which only depends on tq, X, A 
and rZ. Then, for any initial ao satisfying Assumption \2.l\ and coinciding with over 
the boundary dO., the above algorithm generates a seguence {<Jk}, k = 0,1,..., which 
is convergent to and satisfies 

Ikfc - at||L2 < c*^||cro - a^||L2, k = 0,l,..., (5.4) 

where c < 1 depends on ||V(t'^||l°° and hZ. 


Proof We start by subtracting V ■ (crlEfc x Bq) from both sides of (5.2) to obtain 
V ■ ((cTfc+i - at)Efc X Bo) = V ■ (a'f(E'f - E^) x Bq). 

Multiplying the both sides by ak+i — crl and integrating over hZ, we arrive at 


= / ((Tfc+i - cr^)V ■ ((cTfc+i - (T^)Efc X Bo)da; 
Jn 

= [ {ak+i-a^)V-{a^E^-Ek)xBo)dx 
Jn 

= [ (afc+i-at)Va+-((Et-Efc)xBo)dx. 
Jn 
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In the above identities, the first identity follows from a similar argument to the one used 


in the proof of Theorem 4.1 and the last identity follows by noting that (El — E^) x Bq 


is divergence-free. Next, we estimate the electric field difference. As in the proof of 


Theorem 1.1, we conclude that 


|E1 — Efc||j;^2 < C||(jfc — al||2,2. 


By the Cauchy-Schwarz inequality and (5.3), we have that 

< c\\ak 


Ikfc+i 


a' 




and (5.4) follows from an induction argument on k. 


□ 


Remark 5.2 Let us point out that indeed the convergence analysis of the proposed 
algorithm carries through when the inverse problem have a Lipschitz type stability 
In fact, 


estimate. 


Theorem 1.1 still holds true with the condition (1.2) replaced by 
(5.3). Correspondingly, the stability constant will depend on ||VcTl||ioo and blow up as 
IIVXII approaches 2K. 


5.3. Numerical experiments 

Now we present some numerical experiments to verify the convergence theory presented 
in the previous subsection. For each experiment, the true conductivity is assumed 
to be Lipschitz continuous and equal to 0.2 near the boundary and we use constant 
0.2 as the initial model unless otherwise specified. To simplify the computation, we 
transform the 3D problem into a 2D problem by assuming the conductivity is invariant 
along the direction. The setup is as follows. The domain we take is the square 
D = (0,1) X (0,1). We employ a uniform triangulation with a mesh size of 1/64. Both 
the Neumann problem and the stationary advection-diffusion equation are solved using 
a first-order finite element method. The algorithm is implemented using FEniCS, a 
hnite element software package [ 12 ], and using Python as the user interface. All the 
numerical computations are performed on a dual-core laptop computer. 


Example 1. We first consider a simple example. The true conductivity is shown in 
Figure [T^ and the error between the true and reconstructed model is shown in Figure p!bj 
The relative LA-qiioi, —c’'^||/||o'^||, drops to 2.88 x 10“^ after 16 iterations. As shown 
in Figure [T^ a linear convergence rate is observed. 

Example 2. We then attempt to recover a more complicated conductivity model, as 
shown in FigureIn Figure |2b[ we show the absolute error of the recovered model after 
45 iterations. The relative L^-error drops to 2.57 x 10“^ and a linear convergence rate is 
still observed. In this example, the gradient of the conductivity is greater than the one 
in the previous example. According to Theorem [5.1 [ this will lead to a greater prefactor 
c in the convergence rate. The comparison of Figure [T^ and Figure [2c| demonstrates this 
point. 
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Figure 1: (a) true conductivity model, (b) absolute error between the inverted and true 
conductivity model, (c) decay of the relative error in norm (logarithmic scale). 


Example 3. To further investigate and demonstrate the convergence results in 
Theorem 5.1, we perform the third test, which is the “reverse” Example 2. We switch 
the role of the true model and the initial model in Example 2. That is, we try to recover 
the constant conductivity with an initial model as shown in Figure 2a The algorithm 
converges after 1 iteration with the absolute lA-eiioi drops below 5 x 10“®. This implies 
that the prefactor c approaches zero as the true conductivity goes to a constant function. 
Actually, this can be proved by noticing that, when g is constant, the unique solution 


to (5.2) is the same constant for any admissible (Tq. 
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(a) (b) 



(c) 

Figure 2: (a) true conductivity model, (b) absolute error between the inverted and true 
conductivity model, (c) decay of the relative error in norm (logarithmic scale). 


6. Discussion 

We investigated the second step in MAT-MI where the problem is to reconstruct the 
conductivity distribution from internal data obtained in the hrst step. A global Lipschitz 
type stability estimate is established when the conductivity is We devise a 

novel iterative method for solving the inverse problem that involves, at each iteration, 
the solution of a well-posed boundary value problem followed by the solution of an 
advection-diffusion problem. The iterative method is shown to be convergent. Results 
from numerical experiments demonstrate the effectiveness of the approach. 
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It would be interesting to extend the computational method proposed to three 
dimensions and to invert real measured data. An important direction for this research 
is to consider the case of anisotropic conductivity. In [2], the authors examine the effect 
of electrical anisotropy in MAT-MI. A homogeneous tissue is considered. They find 
that, when imaging nerve or muscle, electrical anisotropy has a significant effect on the 
acoustic signal and must be accounted for in order to obtain accurate images. 
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